Emergence of Magnetic Long-range Order in Frustrated Pyrochlore Nd2lr207 with 

Metal-insulator Transition 
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In this study, we performed powder neutron diffraction and inelastic scattering measurements of 
frustrated pyrochlore Nd2lr207, which exhibits a metal-insulator transition at a temperature Tmi of 
33 K. The diffraction measurements revealed that the pyrochlore has an antiferromagnetic long-range 
structure with propagation vector of (0,0,0) and that it grows with decreasing temperature below 
15 K. This structure was analyzed to be of the all-in all-out type, consisting of highly anisotropic 
magnetic moments of magnitude 2.3 ± 0.4^b, where is the Bohr magneton. The inelastic 
scattering measurements revealed that the Kramers ground doublet of Nd^^ splits below Tmi. This 
suggests the appearance of a static internal magnetic field at the Nd sites, which probably originates 
from a magnetic order consisting of Ir*"'' magnetic moments. Here, we discuss a magnetic structure 
model for the Ir order and the relation of the order to the metal-insulator transition in terms of 
frustration. 

PACS numbers: 71.30.-|-h, 75.25.-j, 75.30.-m, 75.50.-y 



I. INTRODUCTION 

Since Pauling's initial proposal for water ice^i geomet- 
rical frustration has played an important role in the fields 
of both chemistry and solid-state physics. Only some 
classical-spin pairs in frustrated magnets can be arranged 
antiferromagnetically on a triangular lattice. ^''^ There- 
fore, frustration suppresses magnetic ordering and pro- 
motes intriguing phenomena in insulators, such as spin- 
liquid-like fluctuations of spin ice and spin molecules,-^ 
and the formation of complex magnetic structures with 
multiferroics j^ii 

Metallic frustrated systems can exhibit novel transport 
properties. Examples of such properties include heavy 
fermion behaviors in spinel LiV204 and in the C15 Laves 
phase Y(Sc)Mn2,^'^ superconductivity in pyrochlore ox- 
ide Cd2Re207,^° anomalous Hall effects in pyrochlores 
Nd2Mo207 and Pr2lr207,llii^ metal-insulator (MI) tran- 
sitions in pyrochlores Cd20s207 and IIg2Ru207, and so 
Q^ ^i3^i4 these systems consist of corner-sharing tetra- 
hedral lattices of magnetic atoms, called pyrochlore lat- 
tices, which provide an ideal platform for the occurrence 
of frustration. Thus, one of the challenging issues in cur- 
rent condensed matter physics is to determine the rela- 
tion between the transport properties and geometrical 
frustration. 

Pyrochlore iridates i?2lr2 07 [R — Nd, Sm, Eu, Gd, 
Tb, Dy, or Ho) also exhibit MI transitions)^""— wherein 
both R and Ir atoms form pyrochlore sublattices. All the 
MI transitions are of the second order . ^^d"^ The transi- 
tion temperature Tmi monotonically increases from 33 to 
141 K with increasing atomic number of the R element 



from Nd to Ho.^^'^^ The R dependence of Tmi suggests 
that the MI transitions are sensitive to the change in the 
ionic radius of R^~^ (i.e. lanthanide contraction) but are 
independent of the 4/ electronic statesJ^iii 

Magnetic susceptibility measurements of i?2lr207 in 
a previous study revealed that the zero-field-cooled 
and field-cooled curves do not follow each other below 
Tmif^*^ This observation suggests that the MI transi- 
tions are accompanied by a magnetic anomaly.^^'^^ Band 
calculations of the iridate in another study revealed that 
the conduction band consists of Ir 5c? and O 2p orbitals 
in the metallic phase^i^ Therefore, the magnetic anomaly 
seems to originate from the Ir magnetic moments. The 
Ir atoms are expected to be tetravalent magnetic ions in 
the insulating phase, as in the case of insulating Sr2lr04 
(Ir''+: 5c^^ Jeff = 1/2, Tcff = 1, low-spin state S = 1/2); 
here, Jeff is the effective total angular momentum, L^s 
is the effective orbital momentum to describe triply de- 
generated t2g states under spin-orbit coupling, S is the 
spin angular momentum, and J^s = icff + S}^ 

Recently, several microscopic experiments have fur- 
ther investigated Nd2lr207, and almost ensured that its 
MI transition originates from magnetic ordering at Tmi. 
For example, Raman scattering experiments of Nd2lr207 
suggested that there is no structural phase transition 
down to 4.2 K and that lattice distortion does not inter- 
estingly cause the MI transition. Further, in studies of 
temperature variation of the Nd 4/-electron states across 
Tmi, crystalline electric field excitations of Nd"^"*" were ob- 
served by inelastic neutron scattering, specific heat mea- 
surements, and magnetic susceptibility measurements, 
and the results were used to deduce the Ir magnetic 
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states JiiSi The ten states ofm^+ (4/^, J = 9/2, L = 6, 
S" = 3/2) split into five Kramers doublets. The first ex- 
cited doublet exists at around 26 meV (300 K), which 
is much higher than Tmi- The ground doublet possesses 
in/out-type Ising moments with 2.37/iB along the (111) 
direction, where fiB is the Bohr magneton. With de- 
creasing temperature, the ground doublet energetically 
splits at Tmi like an order parameter of the MI transi- 
tion, and the splitting energy is about 1.3 meV at 3 K. 
This phenomenon suggests the appearance of a static in- 
ternal magnetic field at the Nd sites, probably caused by 
the magnetic ordering of Ir at Tmi. 

Thus, it is most probable that Nd2lr2 07 simultane- 
ously undergoes the MI transition and the magnetic or- 
dering without structural phase transition. However, the 
following major questions are still open: how the frustra- 
tion is released without structural transition and how the 
MI transition is related to its release. To clarify them, 
information of the magnetic structure, which is prime for 
the study of magnetic frustration, will be necessary. 

In this study, we conducted neutron diffraction and 
inelastic scattering experiments of Nd2lr2 07 in powder 
form. We also studied a magnetic structure model both 
for the Nd and Ir moments below Tmi and determined 
its relation to the MI transition in terms of frustration. 



II. EXPERIMENTAL 

Neutron diffraction experiments were performed on 
the powder difTractometer HERMES (Tl-3) at the In- 
stitute for Materials Research (IMR), Tohoku Univer- 
sity, installed in the thermal guide tube of the JRR-3 
reactor at the Japan Atomic Energy Agency (JAEA)i^ 
Incident neutrons with an initial energy of 24.5 meV 
(A = 1.8204(5) A) were extracted by a (331) reflection of 
a Ge monochromator and horizontal collimation sequence 
of guide-blank-sample-22'. Elastic neutron scattering ex- 
periments were performed on the triple-axis spectrome- 
ter TOPAN (6G) of Tohoku University, also installed in 
the JRR-3 reactor. The final energy of the neutrons was 
fixed to 30.5 meV with a horizontal collimation sequence 
of blank-60'-sample-60'-blank. A sapphire filter and a 
pyrolytic graphite filter efficiently removed fast neutrons 
and higher-order contamination, respectively. 

Inelastic neutron scattering experiments were per- 
formed on the triple-axis spectrometer HER (Cl-1) at 
the Institute for Solid State Physics (ISSP), University 
of Tokyo, installed in the cold guide tube of the JRR- 
3 reactor. The final energy of the neutrons was fixed 
to 3.6 meV with a horizontal collimation sequence of 
guide-blank-sample-radial-blank, where the radial colli- 
mator has three blank channels. A horizontal focusing 
analyzer, which covers the range of scattering angles by 
5° and permits acquisition of higher statistics of data, 
was used. A cooled Be filter and a pyrolytic graphite 
Bragg-refiection filter efficiently removed the half-lambda 
contamination. The energy resolution (full width at half 



maximum) was 0.13 meV under the elastic condition. 

A powder sample of Nd2lr207 was synthesized by a 
solid-state reaction methodJ^iii The sample (4.5 g) was 
placed in a thin aluminum foil and shaped to a hollow 
cylinder with a thickness of 0.7 mm and diameter of 20 
mm in order to mitigate the effect of the strong neutron 
absorption of Ir nuclei as much as possible. Then, the 
cylinder was kept in an aluminum container that was 
placed under a cold head in a ^He or •^He closed-cycle 
refrigerator. 



III. RESULTS 

Figure [TJa) shows the diffraction pattern measured at 
9 K, i.e. below Tmi. In this figure, only fundamental 
reflections of the pyrochlore structure and a tiny peak 
of the Ir02 impurity phase are observed; no superstruc- 
ture lines are detected. Figure [ijb) shows the differ- 
ence between diffraction data measured below and above 
Tmi, i-e. at 9 K and 40 K, respectively. No apprecia- 
ble magnetic diffuse scattering is observed, unlike that 
often seen as the spin-ice state in magnetic pyrochlore 
oxides. ^"^ In contrast, a component characterized by the 
propagation vector Qq = (0, 0, 0) was observed below Tmi. 
Figure [Ijc) shows the (220) fundamental reflection lines 
measured below and above Tmi. A weak yet noticeable 
increase in the intensity is observed as the temperature 
decreases from above Tmi to below it, and no broaden- 
ing of the line width is observed. These results suggest 
a long-range structure of the pyrochlore with Qq. Fur- 
ther, since previous Raman scattering experiments have 
shown that structural transition is not observed down to 
4.2 the <7g structure is most probably magnetic. 

Table U summarizes the difference between integrated 
intensities at 40 K and 9 K for several fundamental re- 
flections in barn/unit cell, where a unit cell refers to a 
crystallographic unit cell denoted by 2(Nd2lr207). The 
scale conversion factor for the units was estimated by 
analyzing the fundamental nuclear reflection intensities 
at 40 K, where the lattice parameter x = 0.33(0) co- 
inciding with those for isomorphic materials Sni2lr207 
and Eu2lr207 as estimated by room-temperature X-ray 
diffraction was obtainedj^i^ The conversion factor also 
includes the absorption correction, since the correction 
is canceled in the ratio to the nuclear reflection inten- 
sities. Further, the reflection- angle dependence of the 
correction can be approximately ignored, since the sam- 
ple shape is cylindrical. We confirmed that the nuclear 
refiection intensities at 40 K are well analyzed without 
the reflection-angle dependence of the correction. 

Figure [Ijd) shows the temperature dependence of the 
summation of the integrated intensities 113 and 222; this 
summation shows the largest intensity increment (see Ta- 
ble With decreasing temperature, the intensity begins 
to increase at around 15 ± 5 K, and not Tmi. A simi- 
lar temperature dependence is observed in the inelastic 
scattering data. Figure [5l^a) shows the energy (E) spec- 
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FIG. 1: (Color online) Measured neutron diffraction data for 
powder Nd2lr207. (a) Diffraction data measured at 9 K. (b) 
Diffraction data obtained by subtracting the 40 K data from 
the 9 K data (a), (c) Bragg reflection lines measured around 
the (220) reciprocal lattice point at 9 K, 40 K, and 102 K. 
(d) Temperature dependence of summation of the integrated 
intensities 113 and 222. In (a), (c), and (d), all the lines are 
a guide to the eye. 



TABLE I: Integrated intensities of magnetic components of 
several fundamental reflections in barn/unit cell. Here, a 
unit cell refers to a crystallographic unit cell denoted by 
2(Nd2lr207). The experimental intensities were evaluated by 
subtracting the integrated intensities at 40 K from those at 9 
K. The best-fit calculated intensities were obtained for the all- 
in all-out structure shown in Figs. [3Jb) and[2I^c). The exper- 
imental and calculated intensities are compared in Fig. |3ja) 
as well. 
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FIG. 2: (Color online) Measured inelastic neutron scattering 
raw data for powder Nd2lr207. (a) Scan data measured at 
constant Q of 0.8 at different temperatureSf2i The ver- 
tical origins shift to the height indicated by the horizontal 
dotted lines, (b) Inelastic scattering intensity distributions in 
{Q,E) space at 3 K. (c) Comparison results of Constant-_E 
scans measured at i5 = 1.2 and 1.4 meV in (b). In (a) and 
(c), the solid curves are a guide to the eye. 



ing behavior is observed with decreasing temperature. A 
broad inelastic signal appears at around Tmi, as shown 
by the dotted arrows; it changes into a sharp peak at 1.1 
meV at around 15 K; below 15 K, the peak sharpens fur- 
ther, as shown by the sohd arrows; in the previous study, 
it was not important for us to know at what tempera- 
ture the peak began sharpening. ■^^ The emergence of the 
sharp inelastic peak can be attributed to the splitting of 
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the ground doublet of Nd 4/ electronsj2i The agreement 
between the data in Figs. [TJd) and [DJa) strongly sug- 
gests that the Qq structure mainly originates from Nd'^'*' 
magnetic moments. 

Figure[5Jb) shows the inelastic scattering intensity dis- 
tributions in (Q, E) space measured below Tui- The hor- 
izontally spreading excitation mode is observed at around 
1.3 meV; it corresponds to the splitting of the Nd'^^ 
ground doublet. The scattering intensity decreases with 
increasing Q, confirming that the excitations are mag- 
netic. However, the mode appears slightly dispersive. To 
verify the dispersion, we compared the constant-£' scans 
measured at 1.2 and 1.4 meV in Fig. [^Jb); the compar- 
ison results are shown in Fig. [2{c). The Q dependences 
at these energies are clearly different, indicating that the 
excitations are not completely flat but dispersive, with a 
bandwidth of - 0.1 meV (1 K). 



IV. ANALYSES 

We analyzed the magnetic structure on the basis of 
the Qo intensities given in Table HI Crystalline field anal- 
ysis2i revealed the magnitude of Nd^+ moments to be 
about 2.37/iB, whereas that of Ir"*"*" moments is expected 
to be IfiB at most. Therefore, as the first approxima- 
tion, we assumed that the Qq intensities consist of only 
the Nd moments. In fact, the statistical errors of our 
data would be too large to resolve the Nd and Ir mo- 
ments. Furthermore, later in this section, the magnitude 
of moment estimated by the magnetic structure analysis 
is confirmed to be in agreement with that estimated by 
the crystalline field analysis. 

The crystalline field analysis also strongly suggests 
that the Nd moments are highly anisotropic along the 
(111) directions (in/out-type Ising moments),— as ex- 
pected from available data on other Nd pyrochlore ox- 
ides. In addition, the measured magnetic susceptibil- 
ity shows a small value of only ^ 10~'^/iB/formula un- 
der a magnetic field of 1 kOe, and no hysteresis curve is 
observed at 5 KJ^^^L^l^ Therefore, among the various 
magnetic structure models described by the in/out- type 
moments, the 3- in 1-out and the 2-in 2-out types are most 
probably ruled out, because these model types are essen- 
tially ferromagnetic. A unique possible solution is the 
all-in all-out type of model, as shown in Fig.[3jb). Thus, 
we examined the consistency between the all-in all-out 
model and the Qq intensities. 

The remaining fitting parameter is the absolute value 
of magnetic moments of Nd'^'*' (m). Therefore, we eval- 
uated it by the least-squares method, by employing the 
magnetic form factor of Nd^+ calculated by Freeman and 
Desclaux.^^ Since the sample was cylindrical in shape, 
the reflection-angle dependence of the absorption factor 
was ignored. The evaluation results showed that the all- 
in all-out model had the best-flt calculated intensities, 
with m(9 K) = 1.3 ± 0.2 /ie (Table |l|; these intensities 
are in agreement with the experimental ones, as shown 



in Fig. Ha). 

Further, from Fig. [H^d), the ratio of magnetic scatter- 
ing intensity at 0.7 K to that at 9 K is estimated to be 
about 3.2. Since the magnetic intensity is proportional 
to the square of m,'^^ the value of m(0.7 K) is estimated 
to be 2.3 ± 0.4 (=^3^ x m(9 K))/iB. This value is in 
good agreement with the value of 2.37/iB estimated by 
the crystalline field analysis>2i, 

To summarize Sees. 3 and 4, we found that the an- 
tiferromagnetic long-range structure with Qq grows with 
decreasing temperature below TncI = 15±5 K. The struc- 
ture can be approximately described by the all-in all- 
out type of model for Nd moments with a magnitude of 
2.3 ± 0.4/iB at 0.7 K. No direct signals for Ir moments 
were obtained in the present experiments. 



V. DISCUSSION 

A. Magnetic structure below Tmi 

The ground doublet of Nd^^^ splits when the temper- 
ature decreases to below Tmi, as shown by the dotted 
arrows in Fig. [2Ka); however, the Nd magnetic structure 
grows mainly below TNd, as shown in Fig. [Ijd) and by 
the solid arrows in Fig. [2ja). This discrepancy in tem- 
perature suggests that not Nd moments exhibit static 
magnetic order at Tmi- In isomorphic Nd2Mo207, the 
magnetic moments of d and / electrons do not order si- 
multaneously, as revealed by neutron diffraction experi- 
ments; the Mo structure grows below Curie temperature 
Tc = 93 K and then the Nd structure grows mainly below 
^ 20 K with decreasing temperaturci^ In analogy with 
Nd2Mo207, Ir magnetic ordering is expected to occur in 
Nd2lr207 at Tmi- The ordered Ir moment is estimated 
to be as small as the errors of m later in this subsection. 

The splitting of the ground doublet of Nd'^^ is slightly 
dispersive, with a width of 0.1 meV (1 K), as shown in 
Fig. Hljb); this width roughly corresponds to the magni- 
tude of Nd-Nd interactions. The isomorphic pyrochlore 
insulator Nd2Sn207 with a similar magnitude of Nd'^'^ 
magnetic moments (2.6/iB) also exhibits a Neel tempera- 
ture of 0.9 K and a Curie- Weiss temperature of —0.3 K, — 
again indicating that the magnitude is about 1 K. Thus, 
the Nd structure growing below TncI is not formed by 
Nd-Nd interactions alone. The structure probably para- 
sitizes the aforementioned hidden Ir order through Nd-Ir 
exchange interactions, whose magnitude is estimated to 
be about the splitting width of 1.3 meV ~ 15 K; this 
width is in excellent agreement with the value of iNd- 

We next discuss the Ir magnetic structure. First, as 
shown in Fig. [2{a), a single sharp peak appears at 3 
K. This fact indicates that the internal magnetic field 
is uniform at all Nd^"*" sites in the entire sample; this 
uniformity in turn suggests that the Ir structure is also 
described by Qq without any modulation. Indeed, Qq has 
no special symmetrical directions and therefore forms no 
domains, and is consistent with the absence of structural 
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FIG. 3: (Color online) Magnetic structure modelling, (a) Comparison between experimental and best-fit calculated magnetic 
reflection intensities given in TableU] (b) All-in all-out magnetic structure. The filled circles represent Nd"^"*" ions, and the arrows 
represent the Nd moments, (c) Trigonally distorted ligand around Ir*+ ion. All the Ir-O bonds are of the same length (2.01 
A), (d) Relation between magnetic moments (blue thin arrows) of Ir''+ ions (blue small balls) and magnetic moments (red thick 
arrows 1 of Nd='+ ions (red large balls) . Both the Ir and the Nd moments form the all-in all-out structures. Alternative directions 
of Nd moments in the case of a ferromagnetic Nd-Ir interaction are depicted: when the moments are antiferromagnetic, the 
directions of all the red arrows are reversed, and the all-in all-out type structure is retained. 



transition. Second, as mentioned above, the antiferro- 
magnetic nature of the magnetic structure is clear from 
the magnetic susceptibility. Third, as shown in Fig. |3^c) , 
the local environment of Ir'^"*' is trigonal, suggesting that 
Ir moments with Joff = 1/2 also favour the (111) direc- 
tions (in/out type). Thus, from these three observations, 
along the same lines as for the Nd structure mentioned in 
§ 4, the Ir moments are also expected to form the all-in 
all-out type of antiferromagnetic structure. 

We can confirm that the all-in all-out type of Ir struc- 
ture generates the all-in all-out type of Nd structure, 
when only nearest-neighbor Nd-Ir exchange interactions 
are taken into account. Figure [21[d) shows the relation of 
Ir and Nd crystallographic sites and magnetic moments 
with the all-in all-out structures. The single Nd site is 
positioned at the centre of the hexagon on the kagome 
plane in the Ir sublattice. At the Nd site, the surrounding 
six nearest-neighbor Ir moments generate the exchange 
field, which is perpendicular to the (111) kagome plane. 
The alternative upward or downward direction depends 
on whether the Nd-Ir interaction is ferromagnetic or an- 
tiferromagnetic. 

To examine the consistency between the approxima- 
tion with only Nd moments in § 4 and the aforemen- 
tioned model with both Nd and Ir structures, it will be 
worth reanalyzing the experimental reflection intensities 
(TableU]), though the statistical errors are large. The fit- 
ting parameters are the absolute value of magnetic mo- 
ments of Nd'^"'" (mNd > 0) and that of (mir) with 
sign corresponding to the alternative upward or down- 
ward direction. We evaluated them by the least-squares 
method, by employing the magnetic form factors of Nd'^"'" 
and f+i^i^ The best-fit calculated intensities (Table IJ) 
were obtained at TONd(9 K) = 1.3 fiB and mir(9 K) = 



—0.1 /iB, which correspond to mNd (0.7 K) = 2.3 /ie and 
TOii(0.7 K) — —0.2 /iB- The negative sign means the 
magnetic structure shown in Fig.[3][c) (ferromagnetic Nd- 
Ir interaction). These calculated intensities are again in 
agreement with the experimental ones within the statis- 
tical errors, as shown in Fig. [31[a). Thus, the very small 
value of mir supports the validity of the approximation 
with only Nd moments. 



B. Relation between frustration and MI transition 

The present experiments revealed that the antiferro- 
magnetic long-range structure with Qq exists below Tmi 
in Nd2lr207. This fact strongly suggests that the frus- 
tration is released below Tmi. That is, the system is ex- 
pected to change from the magnetically frustrated metal- 
lic phase to the antiferromagnetic insulating phase at 
Tmi. The metallic phase will gain orbital hybridization 
energy (band formation energy), whereas the insulating 
phase will gain magnetic ordering energy. 

Then, the following question arises: how is the frus- 
tration released without lattice distortion at Tmi? In 
a pyrochlore lattice, antiferromagnetic nearest-neighbor 
interactions with isotropic moments induce strong frus- 
tration, as is the case in the highly frustrated spinel 
MgCr204 (Cr^+: whereas those with anisotropic 

moments induce the all-in all-out type long-range order 
without frustration. ■^^'■^^ The anisotropy is based on the 
spin-orbit-coupled JeS states. The spin-orbit coupling is 
normally suppressed in a metallic phase (unlike in an in- 
sulating phase), since the coupling constant A is roughly 
proportional to where r is the size of an unpaired 

electron cloud. Thus, in Nd2lr207, above Tmi, an anti- 
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ferromagnetic interaction with relatively isotropic Ir mo- 
ments is considered to hamper magnetic ordering owing 
to frustration. Below Tmi, the spin-orbit coupling will 
become relatively active and the Ir moments will restore 
the in/out-type anisotropy; therefore, the all-in all-out 
type of order will emerge without frustration. 

It is open how the electron/hole bands change below 
and above Tmi. We remark that the Mott-type antifer- 
romagnetic insulator with J^g = 1/2 half-filled bands 
was recently established in Sr2lr04.— The J^s — 1/2 
Mott insulator might be related to the insulating phase 
in Nd2lr207. 

VI. CONCLUSIONS 

In summary, to study a relation among magnetic frus- 
tration, the absence of structural phase transition, and 
the MI transition in Nd2lr2 07, we conducted neutron ex- 
periments of this material in powder form. As the results, 
magnetic Bragg reflections with Qq — (0, 0, 0) were found 
by diffraction. The propagation vector is consistent with 
the absence of structural transition. The reflection inten- 
sity of Nd2lr207 increases with decreasing temperature 
below TNd — 15±5 K. The magnetic structure of this sys- 
tem can be described by the all-in all-out type of model 
for Nd moments with a magnitude of about 2.3 ± 0.4/iB 
at 0.7 K. The magnitude of moments is consistent with 
the results of a previous crystalline field analysis. 

The Kramers ground doublet of Nd'^^ begins to split 
at Tmi with decreasing temperature, suggesting another 
magnetic ordering at this temperature. The magnitude of 



dispersion of the splitting is estimated to be only ~ 1 K, 
which is much below Tnci; this suggests that the Nd mag- 
netic structure cannot be formed by Nd-Nd interactions 
alone. Therefore, we concluded that a hidden magnetic 
order of Ir moments appears at Tmi, and it is parasitized 
by the Nd structure; however, we could not detect direct 
signals for the Ir moments in the experiments, probably 
because of their very small magnitude. 

The Ir magnetic structure can also be described by the 
all-in all-out type of structure. We proposed a scenario 
that the in/out type of anisotropy based on spin-orbit 
coupling plays a key role in releasing the frustration and 
inducing the MI transition. The scenario will provide an 
example of the relation between the transport properties 
and geometrical frustration. 
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